Abstract. Based on SPH simulations, we quantify the 
geometrical distortion effect due to tidal and rotational 
forces on polytropic secondaries in semi-detached bina- 
ries. The main effect is an expansion of the polytropic star, 
with an effect on the radius of ~ 5%-12%, depending on 
the polytropic index and the mass ratio. We apply such 
distortion effects to the secular evolution of secondaries 
in cataclysmic variable systems. We focus on systems be- 
low the 2-3h period gap and that approach the minimum 
period. We find a significant increase of the predicted min- 
imum period (~ 4% if changes in the secondary's thermal 
relaxation are approximately taken into account) . Though 
an improvement, the effect is not big enough to solve the 
mismatch between predicted and observed minimum pe- 
Cj riod at 80 min. 
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1. Introduction 

The description of close binary systems is usually based 
on the Roche model which defines the shape of a binary 
component distorted by tidal and rotational forces. In the 
framework of the Roche model one assumes that the bi- 
nary components (the primary and the secondary) either 
are point masses, or are corotating and have a spherically 
symmetric mass distribution irrespective of their proxim- 
ity or mass ratio (Kopal 1959, 1978). In a semi-detached 
system, one of the components fills its critical equipoten- 
tial lobe defined by the potential of the inner Lagrangian 
point, and which determines the maximum extent of a star 
in a close binary. This is the so-called Roche lobe within 
the Roche model. Cataclysmic variables (hereafter CVs), 
composed of a white dwarf as the primary and a low-mass 
star or a brown dwarf as the secondary, belong to this type 
of systems: the secondary fills its critical lobe and trans- 
fers mass towards the primary. When applying the Roche 
model to problems of binary evolution one makes implic- 
itly the following assumptions (among others). First, the 
Roche potential is a good approximation of the true poten- 
tial that one would obtain by solving the Poisson equation. 
Second, the effects of tidal and rotational forces on the 
internal structure of the star are negligible, i.e. that they 
result in only small corrections compared to stellar models 
assuming spherical geometry. Third, that for the purpose 
of evolutionary computations involving one-dimensional 
stellar models the lobe-filling star may be replaced by a 
spherical star of the same volume. This is tantamount to 
assuming that though tidal and rotational forces change 
the shape of a star they leave its volume invariant. The 
radius of the lobe-filling star then only depends on the 
geometry of the system, and can be calculated by means 
of simple analytical fits (Paczyhski 1971, Egglcton 1983). 
The main purpose of the present paper is to examine in 
some detail the third and to some extent also the second 
of the above assumptions, both of which have so far not 
been tested. 
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Recent 3D simulations (Rezzolla et al. 2001; Motl & 
Frank, priv. comm.) confirm that, at least in the case of a 
semi-detached system, the Roche potential is a good ap- 
proximation if the lobe-filling star is sufficiently centrally 
condensed, i.e. if the effective polytropic index is N > 3/2. 
The analysis of Rezzolla et al. (2001) is based on numeri- 
cal models of semi-detached binaries that account for the 
finite size of the secondary star, thus relaxing the first as- 
sumption inherent in the Roche model. With the validity 
of this approximation for the determination of the poten- 
tial, they also show that such effects hardly affect grav- 
itational quadrupole radiation. Moreover, a comparison 
between the angular momentum loss and mass-transfer 
timcscales predicted by the Roche model and their numer- 
ical models shows small differences. They thus conclude 
that finite size effects cannot account for the mismatch 
between the observed minimum period P m i n at 80 min of 
CV systems and the theoretical value Pturn- The latter is 
indeed ~ 15% shorter than the observed value, according 
to recent calculations based on improved stellar physics 
(see Kolb & Baraffe 1999). Since Rezzolla et al. (2001) 
do not consider thermal relaxation effects in their calcu- 
lations, they can only determine a differential correction 
to Ptum when going from Roche model to self-consistent 
potential. However, in doing so, the second and third as- 
sumptions mentioned above remain untested. The main 
purpose of our paper is thus to explore the consequences 
of making these two assumptions. 

Our main goal is to determine quantitatively the de- 
parture from spherical symmetry of the secondary in semi- 
detached binaries and to analyse the consequences on the 
mass transfer rates and orbital period in CV systems. We 
use smoothed particle hydrodynamics (SPH) techniques 
to study equilibrium configurations of semi-detached bi- 
naries and estimate for different mass ratios the geomet- 
rical deformation of the secondary as it fills its critical 
lobe. The numerical models and results are described in 
§2. In §3, we analyse some of the consequences of the tidal 
and rotational forces on the secular evolution of the low 
mass donor on grounds of models constructed by Kolb & 
Baraffe (1999) and Baraffe & Kolb (2000). We focus on 
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the problem of the minimum period and the discrepancy 
between observations and models (see e.g King 1988 and 
Kolb 2001 for a review on the properties of CV systems). 
A discussion and conclusions follow in §4. 

2. Numerical models of semi-detached binary 
systems 

2.1. The method 

We used a SPH code, originally developed and kindly pro- 
vided by Willy Benz (see details in Benz et al. 1990) to 
perform numerical simulations of a close binary system 
composed of a point mass (primary) and a polytropic star 
(secondary). The SPH method has been described exten- 
sively in the literature (see Monaghan 1992 and references 
therein) and is often applied to the study of close binary 
systems (e.g Benz et al. 1990; Lai et al. 1994; Rasio & 
Shapiro 1995; Segretain et al. 1997) 

In all our simulations, the primary is a 1 Mq point-like 
mass. The secondary is described by a polytropic equation 
of state p = Kp 1+1 / N , where p is the pressure and p the 
density. The polytropic constant K is fixed for a given 
index N by the mass M and radius R of the spherical 
secondary. We adopt two polytropic indices, i.e. N = 3/2, 
which provides a good description of fully convective ob- 
jects such as low mass stars, and N — 3 characteristic of 
solar type stars with M ~ 1M Q . For the particular case of 
CVs, systems below the period gap are well described by 
N — 3/2 polytropes, whereas N — 3 applies to systems 
with periods > 6 h and typical masses around ~ 1 Mq. 
The two values of N thus represent limiting cases for the 
description of CV secondaries. Assuming that K remains 
constant in space applies well to fully convective objects 
with a fully adiabatic structure, and implies a chemically 
homogeneous structure (constant molecular weight) for 
the standard models with N = 3. This is a reasonable 
approximation for the present study. 

The simulations use ~ 15000 particles. In order to 
check the accuracy of our results, we ran a limited num- 
ber of simulations with 57000 particles. We find that 15000 
particles is a good compromise between computational de- 
mand and accuracy. The particles are initially uniformly 
distributed on a hexagonal close-packed lattice. The ini- 
tial number density of particles is constant throughout the 
volume of the sphere describing the initial configuration 
of the secondary. The particle masses are proportional to 
the local mass density. This provides a good spatial res- 
olution near the stellar surface, which is crucial for our 
problem of critical lobe determination where surface ef- 
fects are predominant. The simulations are performed in 
a corotating reference frame with the origin at the center 
of mass of the system. The initial separation A- ln n of the 
two components is arbitrarily fixed at four times the sep- 
aration required for the secondary to fill its Roche-lobe 
^Rocho, estimated from the Eggleton (1983) fit. For such 



a separation, tidal and rotational effects on the secondary 
are negligible. The orbital separation is decreased with the 
arbitrary constant rate (^4i n it — ^4.Roche)/7" s imuj so that the 
total timescale of the simulation r s i mu is ~ 1000 times the 

I r 3 \ 1/2 

typical hydrodynamical relaxation time T ro iax — ( J 
of the secondary. The simulation is stopped when the sec- 
ondary fills its critical lobe i.e. when the first particles 
from the secondary reach the inner saddle point of the po- 
tential []. This marks the onset of mass transfer. Once the 
critical separation is reached, we check that the model has 
reached an equilibrium configuration, starting from such 
critical separation and letting it relax in a non rotating 
reference frame. 

Our goal is to estimate the deformation effects on the 
secondary due to tidal and rotational forces as it fills its 
critical lobe. The deformation can be measured in terms 
of the ratio of the final to initial stellar radius D = Rf/Ri. 
Ri is the radius of the unperturbed spherical polytrope. 
Rf is an effective radius defined as the radius of the sphere 
with the volume Vf of the secondary filling its critical lobe. 
Vf is provided by our SPH simulation at the onset of mass 
transfer. The method to estimate Vf is described in Ap- 
pendix A. 

2.2. Results 

We ran a grid of simulations for various mass ratios 
q = Mi I Mi between the secondary and the primary. Typ- 
ically, CV systems with periods from ~ 10 h down to the 
minimum period cover a range of q between 1 and 0.06. 
Fig. [l] displays the final configuration of a N = 3 polytrope 
with mass ratio q = 0.8. This illustrates the case of CV 
systems with periods > 6 h (see Baraffe & Kolb 2000). 
Fig. H shows the results for the case N — 3/2 and q = 
0.07, characteristic of secondaries approaching the period 
bounce Ptum (Kolb & Baraffe 1999). We note that in the 
case N — 3/2 (Fig. g, lower panel), the surface value of 
<fi is not constant. We did not find any satisfactory expla- 
nation for such behavior. This feature has already been 
noted in some cases by Rasio and Shapiro (1995) and in- 
terpreted in terms of number density of SPH particles be- 
ing not exactly constant around the surface of a star with 
large tidal deformation. Increasing the number of particles 
from 15000 to 57000 and double-checking that the models 
have reached an equilibrium configuration do not solve the 
problem. We do not expect that this affects the accuracy 
of our final results, since our deformation calculations are 
in excellent agreement with similar calculations by other 
authors (see below). 

The resulting deformations D as a function of q are 
summarized in Table 1 for N — 3 and N = 3/2. As ex- 
pected, tidal and rotational distortion yields an expansion 
of the secondary's volume with respect to the unperturbed 
spherical configuration. In terms of effective radius, the ex- 

1 defined as L\ in the Roche potential 
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Table 1. Deformation D = i?f/fi; as a function of mass 
ratio q — M<jonor/-Maccretor of the secondary in a semi- 
detached binary for polytropic indices N = 3 (D 3 ) and 
N = 3/2 (D 3/a ) 



q 


D 3 


D3/2 


0.06 




1.06 


0.07 




1.06 


0.1 




1.05 


0.2 


1.12 


1.05 


0.3 


1.12 


1.05 


0.4 


1.12 


1.04 


0.5 


1.12 


1.04 


0.6 


1.11 


1.04 


0.7 


1.11 


1.04 


0.8 


1.11 


1.04 


0.9 


1.11 


1.04 


1.0 


1.10 


1.04 



pansion is typically 11% for N — 3 and 5% for N = 3/2. 
The dependence of D on N can be understood in terms 
of the compressibility \ = 9 log p/dp — N/(l + N)/p, 
which is larger for N = 3 than for N = 3/2. The larger 
the compressibility, the larger the deformation. 

In order to visualize the deformation of the secondary 
compared to the spherical case, Fig. || displays lines of 
constant density for N = 3 and N = 3/2. An inspection 
of Fig. H shows that the largest departure from spherical 
symmetry is observed in the outermost layers of the poly- 
tropic star, whereas the central regions are only slightly 
affected. The results displayed in Table 1 are in excel- 
lent agreement with the work of Uryu & Eriguchi (1999), 
based on a different numerical method. Indeed, for N — 
3/2 they found distortion effects of ~ 4% for 0.1 < q < 1. 
They however did not analyse the case N — 3. A compar- 
ison of the numerical ratio Rf/A, where A is the orbital 
separation, and the ratio given within the Roche model 
according to Eggleton (1983) shows small differences (less 
than 2%), in agreement with the results of Rezzolla et al. 
(2001) and confirming indeed that the Roche potential is 
a good approximation in the present case (the so-called 
first assumption, see §1). 

3. Application to cataclysmic variable 
evolutionary models 

In order to analyse the consequences of the distortion ef- 
fects found in the previous section on period and mass 
transfer rate in CV systems, we follow the secular evolu- 
tion of the secondary using the same models and input 
physics as described in Kolb & Baraffe (1999) and Baraffe 
& Kolb (2000). We focus on systems below the 2-3 h period 
gap and specifically on the minimum period discrepancy 
between observations and models. Although distortion ef- 
fects seem to be more important for systems above the 
period gap (P > 3 h) (see Tab. [l]), their consequences 



are difficult to quantify given the large uncertainties of 
evolutionary models describing such systems, such as the 
magnetic braking law and the resulting mass transfer rate, 
the evolutionary stage of the secondary at onset of mass 
transfer or the mixing length parameter. Below the pe- 
riod gap, such uncertainties are fortunately considerably 
reduced (see Baraffe & Kolb 2000; Kolb et al. 2001 for 
details) . 

3.1. Geometrical effects 

In the following, we only consider the effects of distor- 
tion on the geometry of the system. The orbital proper- 
ties, e.g. the orbital period and separation, and the mass 
transfer rate will be indeed affected by the larger effec- 
tive radius of the donor, estimated in §2, compared to the 
undistorted case. However, for the moment, we ignore the 
rotational and tidal effects on the thermal structure of the 
star, assuming that its inner structure is unaffected and 
determined by the unperturbed stellar structure equations 
in spherical symmetry. A rough estimate of the thermal 
effects on the secondary's properties resulting from its ex- 
pansion is derived in the next section (§3.2). 

We analyse an evolutionary sequence with an initial 
donor mass of 0.21 M®, a primary mass of 0.6 M@, and 
gravitational radiation (GR) as angular momentum loss 
mechanism (see Kolb & Baraffe 1999). From the radius 
i?2 obtained from integration of the standard stellar struc- 
ture equations, and the mass ratio q, which varies along 
the sequence of evolution, the effective radius is derived 
according to Tab. 0. The mass transfer rate is then cal- 
culated as a function of the difference between effective 
donor radius and Roche radius, following Ritter (1988). 

The comparison between sequences without distortion 
(solid line) and with distortion (dashed line) is shown in 
the orbital period - effective temperature diagram (Fig. 
[|). Although reducing the discrepancy with the observed 
Pmin, distortion effects provide an increase of the mini- 
mum period Ptum of only ~ 6% (or ~ 4-5 min), com- 
pared to the undistorted case. This is slightly less than 
what is naively expected from the period - radius rela- 
tion P oc (R^jMi) 1 ^ 2 . An increase of the radius by ~ 
6%, as expected from distortion effects near the minimum 
period (see Tab. Q), should indeed yield ~ 9% increase 
of P. The smaller effect found on P stems from the de- 
pendence of angular momentum loss driven by GR on 
the secondary radius Jgr/J oc P~ 8 / 3 oc R^ 4 - Conse- 
quently, the larger radius in the distorted sequence im- 
plies a decrease of JgRj and thus a smaller mass transfer 
rate —Mi. As shown below Ptum depends on the ratio 
t = tKn/tM of the secondary's Kelvin-Helmholtz time 
£kh and the mass transfer time tu — Mi/ (—Mi), The 
decrease of —Mi thus yields a decrease of r, implying less 
departure from thermal equilibrium and thus a smaller 
-Ptum- Because Jqr depends explicitly on the mass of the 
primary, Ptum does also depend on it, but only weakly, as 
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Fig. 1. Configuration of a polytrope N = 3 and mass ratio q — 0.8 as it fills its critical lobe, i.e at the onset of particle 
transfer towards the primary. The upper panel is a projection of the SPH particles onto the orbital x — y plane, z 
being the rotation axis, x is the coordinate along the binary axis. The primary is indicated on the left side of the 
plot by a thick dot. The lower panel displays the projection onto the (x, <£>) plane, where $ is the effective potential 
i.e the sum of the gravitational and the centrifugal potentials (rotating reference frame). The solid line indicates the 
variation of Q(x,y = 0, z = 0) along the binary axis. In both panels, the inner Lagrangian point is indicated by a 
symbol surrounded by an open circle. The simulation is done with 15000 particles. 



shown by Paczyhski & Sienkicwicz (1983), who found that 
<91nP t urn/<91n Mi « 0.09. In fact our computations show 
that Ptum varies from 71 min for Mi = 0.6 M Q to 74 min 
for Mi = 1.2 Mq, when distortion effects are included. 

3.2. Thermal relaxation 

In the numerical computations discussed in the previous 
section, we have not taken into account changes in the 



thermal reaction of the secondary which must result from 
its inflation due to tidal and rotational forces. A fully 
consistent treatment of the distortion effects would imply 
solving the multi-dimensional stellar structure equations. 
Rather than doing this we shall in the following derive a 
rough estimate of the thermal relaxation effects and ex- 
plore their consequences for the minimum period of CV 
systems. We indeed expect that the changed surface area 
of the more distended secondary, as a result of the distor- 
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Fig. 2. Same as Fig. [j] for a N = 3/2 polytrope and mass ratio q = 0.07. The simulation is done with 57000 particles. 



sion effects, will affect its surface luminosity, and thus its 
thermal properties. 

We denote by i?2,o, £2,0, T c s,q, ■ ■ ■ the quantities of the 
donor star which result from assuming a pure 1/r poten- 
tial, and by i?2, £2, Tcff, • ■ • the corresponding quantities 
of the spherical equivalent of the critical lobe-filling star. 
Obviously we have 

i?2 = R.2,0 -Dn , (1) 
where is the deformation factor for a polytropic index 
N determined in §2.2 (see Table 1). We thus have for the 
orbital separation A and the orbital period P 

A = Aq -Dn , (2) 
and 



P = PoD. 



3/2 



(3) 



Since mass transfer below the period gap is assumed to 
be driven by loss of angular momentum via gravitational 
radiation alone, the mass transfer rate can be written as 
(see e.g. Ritter 1996) 



Mo = 



M 2 



-2 



Jgr 



(4) 



Ccff - CcL \ J J 

Here Coff = d In R2 /d In M2 is the effective mass radius 
exponent of the donor star, and £cl the mass radius ex- 
ponent of the volume-equivalent critical lobe radius^. Be- 
cause 

>^gr -4 



J 



oc a 



(5) 



2 Ccl is equivalent to the mass radius exponent of the Roche 
radius £r in the Roche model (see e.g. Ritter 1996) 
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Fig. 3. Lines of constant density (in log of the secondary 
in the x — y plane: dotted lines correspond to the spheri- 
cal unperturbed secondary and solid lines to the distorted 
case. Upper panel: N = 3 and g=0.8. Lower panel: N = 
3/2 and g=0.07. The mass centers of the spherical and 
distorted polytropes are located at (x,y) = (0,0). 



we have 

-Ma = -A>2,o£n~ 4 , 
assuming that £ G g- is the same. 
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Fig. 4. Effective temperature versus orbital period for evo- 
lutionary tracks with undistorted (solid line) and distorted 
(dashed line) secondary. The dash-dotted line corresponds 
to the case with distortion, including thermal relaxation 
effects as estimated in §3.2. A test case with distortion 
effect included and J — 2.5 Jgr is also shown (dotted 
line). 



Because of the larger radius of the critical lobe-filling 
star, its surface is larger, thereby affecting its luminosity, 
effective temperature and Kelvin-Helmholtz time 



^kh 



GMo 



GM 2 



47raR 2 3 T eS 4 ' 



(7) 



R2L2 

Consequently, its thermal relaxation, i.e. its gravo- 
thermal luminosity L g will also be different. 

In order to estimate the change of ^kh, we need to 
determine the change in T e g with radius i?2- Since donor 
stars below the period gap are fully convective, we can 
apply the theory of the Hayashi-line, as described by Kip- 
penhahn and Weigert (1990). Accordingly we get 

log R 2 



l0gT cff = 



5a 



-26- 
a + 



logM2 -I- const. 



(8) 



5a + 26 + 5 

Here a = (d\ogK,/d\ogp) T and 6 = (d\ogn/d\ogT) p , 
where k is the photospheric opacity. Typically, for very 
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low-mass stars with atmospheric opacities dominated by 
molecular absorption, we have a rj 1 and 6w0. Therefore 

aiogTeff^ _ 1 



(9) 



M 



d log R 2 

We can now use ([[]) together with (JsJ) in (0) and obtain 

^KH = ^KH.O Dj$ ' . (10) 

We note that in deriving (|o|) we have applied in (Q) the 
factor _Dn only to the P-dependence coming from the lu- 
minosity Indeed, the R 2 in the denominator of (0) 
comes from the gravitational binding energy of the star 
which is related via the Virial theorem to its thermal en- 
ergy, i.e. essentially to its central temperature. Since the 
central regions of the star are expected to be hardly af- 
fected by the distortion of the outer layers (cf. Fig. 3), it 
is not appropriate to also propagate the factor P>n to the 
remaining factor R 2 . 

Let us now examine the conditions at the period Pturn- 
For a polytope of index N losing mass the effective mass 
radius exponent can be written as (e.g. Ritter 1996) 
r _ > , 5 ~ N *m 2 L e 

Cefi — Cad + 7, ~ 7 — ' 



where 



3 - N i KH L 



c 



ad 



dlogP 



1 - N 



(12) 



<91ogMy K 3 — N 
is the adiabatic mass radius exponent and 

f M 2 = r*~ (13) 
M 2 

is the mass loss time scale. Because at P = Pt urn the 
donor star is characterized by TV = 3/2 and £ c ff = +1/30 
we obtain from ( pi] ) and ( fL2| ) 
*M 2 _ 2 
£kh L 7 
at P = Pturn, or 



2 ^KH 
7 ^M, 



2 tKH 
7 ^M 2 



3/2 



-6.8 



(15) 



'3/2 

0.67 



(16) 



With 73.3/9 w 1.06 (cf. Table 1) (gj) yields 
L V L 

Eq. ( |l6| ) means that the 3D-effects, i.e. the reduced mass 
loss rate and the increased surface luminosity of the more 
distended star, result in a smaller deviation from thermal 
equilibrium at P = Pt U m ■ Hence the star is systematically 
less inflated by the effects of thermal disequilibrium, and 
this, in turn, compensates at least partially for the sys- 
tematic increase of the orbital period due to the factor 
D a /2 > 1. 

A quantitative estimate of the decrease of Pturn sug- 
gested from (^) can be derived by recomputing the evo- 
lutionary sequences including the effect of distortion, as 
done in §3.1, and by artificially increasing the radiating 
surface of the donor in the Stefan-Boltzmann law by a 
factor -D3/2 2 ' 8 , as suggested from (10). Note that this is 
equivalent to increasing the radiating surface by a factor 



3 this follows directly from Kepler's 3rd law, the fact that 
(1 + l/q)V f /A 3 w const, for g<0.8, and P = 



-D 3 / 2 2 , and to reducing the surface gravity in the inte- 
gration of the stellar atmosphere by the same factor. The 
result of such a numerical experiment is displayed in Fig. 4 
(dash-dotted line) and shows a slight decrease of Pt U m by 
~ 2-3% compared to the case with pure geometrical effects 
(dashed line). These results fully confirm the expectation 
derived from (16), namely that the value of Pturn is re- 
duced by taking into account the effects of the changed 
thermal relaxation. 



4. Discussion and conclusions 

Although reducing the discrepancy between observed and 
predicted minimum period, distortion effects seem insuf- 
ficient to provide a satisfactory solution of the mismatch 
between calculated and observed minimum period. A com- 
bination of distortion effects as estimated in §2 and an 
angular momentum loss rate of 2-2.5 x Jqr can reconcile 
Pturn with the observed 80 min value (see dotted line in 
Fig. 0). Note that without distortion effects, one would 
need 4 x Jqr to reach ~ 80 min, as estimated in Kolb & 
Baraffe (1999). The more modest increase of J required 
according to our calculations is also in better agreement 
with Patterson's (1998) estimate based on space density 
considerations. Additional physical processes can also re- 
sult in an inflation of the secondary, e.g. irradiation from 
the primary (Ritter et al. 2000) or star spots (Spruit & 
Ritter 1983). A rough estimate of irradiation effects or 
star spots can be derived by following Ritter et al. (2000), 
i.e. by reducing the effective radiating surface of the star 
by a factor (1 — s c g). Adopting in our secular evolution 
calculation a factor s c s — 1/2 and J = Jgr yields a se- 
quence very similar to the one obtained with deformation 
(dashed line in Fig. ||). If s c ff = 2/3 the result resem- 
bles the sequence with deformation and J = 2.5 x Jq R 
(dotted line in Fig. 0) . In order to know whether such val- 
ues of Scff, yielding effects comparable to the distortion 
effects, are reasonable requires a sophisticated treatment 
of star spots or irradiation. An investigation of irradia- 
tion effects on non-gray stellar atmospheres is in progress 
(Barman 2001). We stress however that even if distortion, 
irradiation, star spots or additional sources of J are possi- 
ble solutions for removing the mismatch between observed 
and predicted minimum period, the so-called period spike 
problem still remains. A period spike which is a conse- 
quence of the accumulation of systems near Pt U m (where 
P = 0) is indeed predicted by all models for which J or 
Soft are assumed to be the same for all systems. Even if 
they are not, it is very difficult to "smear out" the period 
spike in a population of systems with different individual 
bounce periods (Barker & Kolb, in preparation). Such a 
period spike is, however, not observed (see Kolb & Baraffe 
1999). 

Finally, we note that Kolb & Baraffe (1999) obtained 
negligible effects on the secondary structure and evolution 
when applying tidal and rotational corrections to the ID 
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stellar structure equations, on the basis of the scheme by 
Chan & Chau (1979). The effect on the total radius in 
Kolb & Baraffe (1999) is much smaller (< 2%) than that 
found from the present SPH simulations. Since two dif- 
ferent numerical methods, on the one hand the work by 
Uryu & Eriguchi (1999) and on the other hand the present 
work, predict the same quantitative deformation effects, 
we are confident that the results of our SPH simulation 
are accurate. Although in the SPH simulations we do not 
take into account the thermal reaction of the star to its 
inflation and deal with polytropes, on the basis of our sim- 
ple estimate given in §3.2 we do not expect the thermal 



effects to significantly reduce the radius of the deformed 



star. A possible reason for the discrepancy between the 
calculations by Chan & Chau (1979) and the present re- 
sults could be the limitation of the former ID scheme to 
describe multi-dimensional effects. Figure 3 shows strong 
effects in the outermost layers which may be difficult to 
account for with such a scheme. In any case, both ap- 
proaches have their shortcomings, but they both provide 
the same conclusion regarding the mismatch of the ob- 
served and predicted minimum period. 

To conclude, our SPH simulations suggest that tidal 
and rotational distortion effects on the secondary in semi- 
detached binaries may not be negligible, and may reach 
observable levels of ~ 10% on the radius for specific cases 
of polytropic index and mass ratio. Although this effect 
yields an increase of the predicted minimum period for 
CV systems, it remains too small to explain the observed 
value of 80 min. Additional effects such as irradiation, star 
spots or extra sources of angular momentum loss still seem 
to be required, leaving the problem of the minimum period 
of CV systems unsettled. 
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Appendix A: Determination of the volume in the 
SPH calculation 

To estimate the final volume of the secondary at the onset 
of mass transfer in our SPH simulations, we proceed as 
follows. We first determine the smallest rectangular box 
containing the secondary star. This box with volume Vbox 
is then filled with iV to t points following a Sobol sequence of 
pseudo-random numbers. Such a sequence is self-avoiding, 
i.e. the points are spread out randomly but in a uniform 
way (see Press et al. 1992), allowing a more uniform filling 
of a volume than a standard random method. Each point 
within twice the smoothing length of a particle is counted, 
otherwise it is discarded. The number of points iV ffective 
fulfilling such condition determines the volume V^tar of the 
star: 
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1 



A tot 

The radius of the star i? s t a r is then defined as the ra- 
dius of the sphere of same volume: 



This method has been tested on spherical and elliptical 
structures of known volumes. According to these tests, it 
provides the radius within the size of one particle i? pa rt- 
Typically, for N part — 15000, we have 



R 



part 



N, 



1/3 
part 



24 



Our method thus determines the radius of a star within 
a systematic error of 4%. Although rough for a precise 
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determination of a stellar radius, the uncertainty is much 
smaller on the ratio of the radius of the same star at two 
different times of the simulation. This is the case for the 
deformation D which is the main quantity of interest in 
our analysis. If e is the absolute error on the radius, R\ 
and i?2 the stellar radii at respectively time t\ and £2, one 
can write: 



R2 + e R 2 I \ ( 1 1 



i?i + e R\ V V R2 Ri 



+ o(e) 



For the specific case of distortion calculation, i?i is the 
radius of the unperturbed polytrope and R2 the radius of 
the polytrope filling its critical lobe. We find typical val- 
ues of R2/R1 < 12% (see Tables Q) with e ~ 4%. Thus 
the first-order term in our last equation provides a correc- 
tion of at most 0.5%. Consequently, our method used for 
estimating the volume is accurate enough for the present 
study. Note that the test simulations done with 57000 par- 
ticles give the same distortion factor (by less than 1%) 
than calculations with 15000 particles. 



